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SUMMARY 


This report provides the analytical documentation for the SOURCE3D Rotor Wake/Stator 
Interaction Code. It derives the equations for the rotor scattering coefficients and stator source 
vector and scattering coefficients that are needed for use in the TFaNS (Theoretical Fan Noise 
Design/Prediction System). Source vector coefficients are the coefficients for stator scattering of 
rotor wakes into acoustic and vortical waves. Scattering coefficients are the reflection and 
tr ansmis sion coefficients for the rotor and stator for inputs of unit amplitude acoustic or vorticity 
waves. SOURCE3D treats the rotor and stator as isolated source elements. TFaNS uses this 
information, along with scattering coefficients for inlet and exit elements, and provides complete 
noise solutions for turbofan engines. 

In the past, three-dimensional source codes, generally, have only treated noise produced 
by the rotor wake interacting with the stator. Reflection and tr ansmis sion from the rotor has been 
ignored. Also, vorticity waves and the effect of swirl have been omitted. These items are needed 
to account for frequency scattering and mode trapping. Frequency scattering is the phenomenon 
where the rotor scatters an input wave at blade passing frequency (BPF) into other harmonics. 
Mode trapping is the mechanism whereby the BPF mode reflects back and forth between the rotor 
and stator, then finally scatters into higher harmonics, and exits. 

SOURCE3D is composed of a collection of FORTRAN programs that have been obtained 
by extending the approach of the earlier V072 Rotor Wake/Stator Interaction Code. Similar to 
V072, it treats the rotor and stator as a collection of blades and vanes having zero thickness and 
camber contained in an infinite hard-walled annular duct Acoustic and vorticity waves are 
calculated by distributing pressure dipoles on the blades and vanes. Then acoustic wave 
amplitudes are calculated at arbitrary points within the duct using the normal mode expansion of 
the Green’s function for an annular duct; vorticity wave amplitudes are calculated via a strip 
application of S. N. Smith’s cascade theory. 

SOURCE3D adds important features to the V072 capability - a rotor element, swirl flow 
and vorticity waves, actuator disks for flow turning, and combined rotor/actuator disk and 
stator/actuator disk elements. These items allow reflections from the rotor, frequency scattering, 
and mode trapping, thus providing more complete noise predictions than previously. The code 
has been thoroughly verified through comparison with D. B. Hanson’s CUP2D two- dim ensional 
coupled cascade noise prediction code using a narrow annulus test case. For this case, which 
provides a two-dimensional situation, our results were essentially identical. 
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CHAPTER 1 


INTRODUCTION 

In turbofan noise prediction programs such as the V072 Rotor Wake/Stator Interaction 
Code (Ref. 1), rotor wakes impinge on the stator producing unsteady vane loads that cause noise. 
Generally in these codes, reflections from other fan elements are ignored. However, reflections 
from the inlet, exit, and rotor could affect source strength. Experimental studies by Topol, 
Holhubner, and Mathews (Ref. 2) have shown this to be the case. Also, Hanson has 
demonstrated this analytically using his coupled cascade theory (Ref. 3) and Meyer (Ref. 4) has 
shown that inlet reflections coupled with the stator can have an impact on source noise levels. It 
was shown in Ref. 3 that such such reflections are necessary for accurate noise predictions; also 
that a swirl component of mean flow and vorticity waves must be included. These additions to 
the acoustic models account for mode trapping and frequency scattering which contribute to 
improved correlation with empirical results. Mode trapping is the phenomenon, described by 
Hanson (Ref. 3), where a mode at blade passing frequency (BPF) reflects back and forth between 
the rotor and stator, but is cut off in the inlet and exit Scattering, i.e. frequency scattering, of this 
mode by the rotor produces higher harmonic modes which then escape through through the inlet 
or exit The presence of swirl between the rotor and stator provides a speed range where mode 
trapping occurs. 

The TFaNS (Theoretical Fan Noise Design/Prediction System, Ref. 5), developed by Pratt 
& Whitney under contract to NASA Lewis, seeks to overcome the shortcomings of earlier noise 
prediction systems by including unsteady coupling between four acoustic elements - inlet rotor, 
stator, exit - shown in Fig. 1.1. It extends the two-dimensional approach developed by Hanson in 
Ref. 3, coupling the fan elements via the duct eigenmodes. 

In TFaNS, the acoustic elements are coupled at the three interface planes shown in 
Fig. 1.1. The axial coordinate is x } and the axial locations of the interface planes are at x 1 , x 2 , x 3 . 
At these planes, state vector A with elements A 1 , A 2 , and A 3 , and source vector B with elements 
B 1 , B 2 , and B 3 , are defined. The state vector elements A 1 , A 2 , and A 3 give the pressure modal 
amplitudes for acoustic waves and velocity modal amplitudes for vorticity waves. The source 
vector elements B 1 , B 2 , and B 3 give the pressure modal amplitudes for the acoustic waves 
generated by the rotor wake impinging on the stator. These elements also give the velocity modal 
amplitudes for the vorticity waves generated by this wake. State vectors and source vectors will 
be described later, in more detail, in Section 2.2. 

The action, in TFaNS, of the rotor wake impinging on the stator to create waves that pass 
back and forth between the various elements to produce the final noise, is described by the system 
of equations 


A = SA + B, 


( 1 . 1 ) 


- 1 - 




1 

2 

3 

INTERFACES PLANES 

X 1 

x2 

x3 

AXIAL COORDINATES 

A 1 

A2 

A3 

STATE VECTOR ELEMENTS 

B 1 

B2 

B 3 

SOURCE VECTOR ELEMENTS 


FIGURE 1.1 ACOUSTIC ELEMENTS FOR A TURBOFAN ENGINE 



where S is a scattering matrix made up of modal reflection and transmission coefficients at the 
various interfaces. The scattering coefficients will be described in more detail later. Solving this 
equation gives the state vector 




(1.2) 


which can then be used to calculate sound power levels and other desired acoustic information. 

At the heart of the TFaNS system is the SOURCE3D Rotor Wake/Stator Interaction 
Code which calculates the matrix S and vector B needed for the equations above by treating the 
stator and rotor source elements each in isolation. The purpose of this report is to document the 
approach used by SOURCE3D in attaining this goal. 

SOURCE3D was developed at Hamilton Standard; it is coded in FORTRAN 77 and was 
designed for UNIX™ workstations. It uses the basic framework provided by V072 and retains 
much of that code. The basic physical and geometric assumptions for SOURCE3D are the same 
as those for V072. In particular, SOURCE3D treats unsteady, subsonic, isentropic flow in an 
annular, infinit e duct There is a background of uniform mean flow in each of the element 
regions. The stator vanes and rotor blades are treated as flat plates with negligible camber and 
thickness. Ref. 1 provides further details regarding basic assumptions. 

SOURCE3D is significantly different from V072 in that it treats both the stator and the 
rotor, rather than just the stator. Also, it calculates source vector coefficients and scattering 
coefficients, rather than the pressure modal amplitudes computed by V072. Additionally, solid 
body swirl has been added to the mean flow between the rotor and stator, along with actuator 
disks to implement the flow turning accompanying the swirl. These disks turn the flow at the 
upstream edge of the rotor and straighten it at the furthest trailing point of the stator. Only 
acoustic waves were included in V072. SOURCE3D adds vortidty waves and the source vector 
coefficients and scattering coefficients associated with these. Rather than treat the stator and its 
actuator disk as separate acoustical elements, SOURCE3D first treats these individually, then 
combines them into a single element The same is true for the rotor and its actuator disk. 
Acoustic wave source vector coefficients and scattering coefficients are computed in much the 
same way that V072 (Ref. 1) determines modal pressure amplitudes. It calculates vorticity wave 
type source vector coefficients and scattering coefficients much the same as modal amplitudes are 
computed in Smith (Ref. 6). In order to provide the features mentioned above, much new code 
and numerous new subroutines have been added to V072. Additionally, many portions of the old 
code have been modified. 
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CHAPTER 2 


PRELIMINARIES 


Before deriving, in the next chapter, source vector coefficients and scattering coefficients 
for the stator and rotor, we need to provide some preliminary information. The first section of 
this chapter treats geometry and provides a definition of swirl by means of a coordinate system 
transformation. Much of this is repeated from Ref. 1, but in less detail; the swirl part is new. 
Standard waves, source vector coefficients, and scattering coefficients are defined in the second 
section. The standard waves provide the form of the input and output for scattering. They are 
also treated in Ref. 5. There is additional information about source vector coefficients and 
scattering coefficients in Refs. 3 and 5. 


2.1 Geometry and Definition of Swirl 

The rotor/stator geometry is shown in Fig. 2.1. As seen there, a cylindrical polar 
coordinate system is fixed in the duct with axial coordinate xj along the duct centerline. The 
coordinate xi increases in the direction of air flow. Depending on the context, the origin for xj 
may be different in different situations. The rotor turns with fan rotational speed Q in the 
direction of increasing <|>. The duct has inner radius r s and outer radius r D . The rotor is made up 
of B evenly spaced identical blades, and the stator consists of V evenly spaced identical vanes. 

The axial and azimuthal sweep of the rotor blades and stator vanes are defined, 
respectively, by the parameters xrled, ymxD, and xsd, yso specified in Fig. 2.2. These quantities are 
defined relative to the leading edges of the rotor and stator at the hub. For the definition of the 
wake, the parameter yso, shown in Fig. 2.2 and defined relative to the rotor tr ailin g edge at the 
hub, is needed. The parameter xspac, also required for the wake, corresponds to the axial spacing 
between the rotor trailing edge and the stator leading edge. All these parameters are functions of 
radius r and are zero at the hub. 


Fig. 2.3 shows an “unrolled” rotor/stator geometry, i.e., the intersection of the rotor 
blades and stator vanes with a cylindrical surface of radius r. In this figure, as is the stagger angle 
of the stator and a* is the negative of the stagger angle of the rotor. Angle a R is defined this way 
so that it will be the analog of as when switching from stator- to rotor-based derivations. The 
blades have spacing 2nr/B, and the vanes, InriV. The local semi-chords of the blades and vanes 
are, respectively, b R and b which are functions of r. Note that the (xj, x 2 )-a xes, shown fixed to the 
stator, will be used later for analysis with regard to the stator. A similar coordinate system, 
shown fixed to the rotor, will be used later for rotor analysis. Other parameters included on the 
figure are there for later reference and will be introduced then. Also note that some of the 
parameters on Fig. 2.3 are further defined on Fig. 2.2. 
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FIGURE 2.2 BLADE/VANE GEOMETRY 



ROTOR L.E. 
RADIUS=r 



FIGURE 2.3 ROTOR/STATOR SCHEMATIC 
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For the case of acoustic waves, solid body swirl is approximated by rotating the stator into 
the approaching flow. It is introduced by using the same three-dimensional approach as in Ref. 1 
only allowing both stator and rotor to rotate with speeds, respectively, of Q; and (see 
Fig. 2.3). The difference in speeds corresponds to the actual rotor speed Q. The stator speed Q.i 
corresponds to -Q.„ where £2, is the swirl speed.* Swirl is input by evaluating source vectors for 
the stator and scattering coefficients for both the stator and rotor in the stator-fixed coordinate 
system rather than in the inertial one. For the case of vorticity waves, the use of different 
coordinate systems is not necessary, because swirl is included direclly, using a two-dimensional 
approach similar to that of Smith (Ref.6). 


22 Standard Waves and Source Vector/Scattering Coefficient Definition 

Three types of waves exist for two-dimensional duct flow (Refs. 3 and 6), upstream- and 
downstream-going pressure waves, and downstream-going vorticity waves. For three- 
dimensional flow, the physics is much more complicated and is the subject of current research 
(Ref. 7). However, here we take the wave families to be the upstream- and downstream-going 
pressure waves and two independent vorticity waves convected with the mean flow. The first 
vorticity wave can be defined as a wave with axial and tangential velocity components (no radial 
component) and corresponds to the two-dimensional vorticity wave; the second can be defined as 
having tangential and radial components (no axial component). The pressure waves propagate at 
the speed of sound; the vorticity waves are convected at the mean flow speed and do not create a 
pressure disturbance. At this stage of development, SOURCE3D uses only the two acoustic 
waves and the first vorticity wave. 

These three types of waves are represented in terms of standard wave forms, which in turn 
are specified in terms of modal amplitudes, or applying more general te rmin ology, state vectors. 
Using these quantities, we can define source vector coefficients and scattering coefficients. The 
standard wave forms presented in this section are those that pertain to the stator-fixed coordinate 
system which is the system in which all our final results are given. 

Pressure waves, either upstream- or downstream-going, are represented by the expression 


111 

5=— oo fc= — oo n=l 


A P P , s I - 7 ( X] - j : P ) -sBQl 1 
A WsknVmn( r )e 1 J > 


( 2 . 1 ) 


where the state vector components, A^ shl , represent the pressure modal amplitudes and are 
normalized by the far-field pressure, /?« . This expression holds in the vicinity of the axial 
interface locations x p , where P = 1, 2, 3 (see Fig. 1.1). The subscript W = 1 represents 
upstream-going pressure waves, W=2 downstream-going ones. Further, s is the blade passing 

* The reason for the minus is that swirl is defined as positive in the direction of positive x 2 
(Fig. 2.3). However, retaining previous V072 conventions, <|> and Q; are positive in the direction 
of negative x 2 . 
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frequency (BPF) harmonic index; k, the vane passing frequency harmonic index; n, the duct 
radial mode index;* and m is the circumferential mode order, where m = sB - kV. The function 
Y £ n (r) is the duct radial mode, y wskn is the axial wavenumber,** and t is time. In the code, as 

would be expected, the infini te limits for s, k, n, as seen in Eq. (2.1), are truncated to give a finite 
number of propagating and decaying modes for the coupling analysis. The user sets the criteria 
for how this is done. 

Note that not all of the indices appearing in the Ref. 5 version of Eq. (2.1) are present in 
Eq. (2.1), because the same degree of generality is not required here. Further, note that, as we 
proceed through this report, for convenience, we may sometimes omit some indices from symbols 
such as those in Eq. (2.1), when it is clear from context what these indices should be. Also, the 
variable x P , when indicated, may sometimes represent locations other than those in Fig. 1.1. 

The pressure in Eq. (2.1) is a linear combination of the standard modes for an annular duct 
with hard walls. These modes are discussed in Ref. 1 and have the general form 




( 2 . 2 ) 


The radial modes here and in Eq. (2.1) are given by 

Y raiW = Arm^m( K jnn r ) » (2-3) 


where J m and Y m are Bessel functions of the first and second kind, and K™ is the duct mode radial 
eigenvalue. The duct mode amplitudes A™ and Bm, are chosen so that the ma xim u m value of Y™ 
over r is +1. Further, the axial wave numbers are given by 


1 



( sBQ. - rnQ.. 'I 

.. 1 

M 

l+ 

f 

V c o ) 

J 


(2.4) 


where the (+)-sign is for upstream-going waves (W = 1) and the (-)-sign is for downstream-going 
ones (W = 2). Also, c 0 is the speed of sound; M is the axial flow Mach number, Ulc 0 , where U is 

the duct unif orm axial velocity; p = Vl - M 2 ; and 


* Consistent with the convention in V072, radial mode indices here start at 1. However, in 
conformance with standard practice, when output is printed, SOURCE3D adjusts this index to 
begin instead at 0. 

** These wavenumbers are the negative of those specified for TFaNS in Ref. 5. However, the 
source vector coefficients and scattering coefficients calculated will be the same because the sign 
in front of y^fot ® Eq. (2.1) is minus (-) and it is plus (+) in TFaNS documentation. 
SOURCE3D utilizes the present convention in order to maintain consistency with V072. 
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(2.5) 


^ skn 


r 

sBQ. - m£l s 

co 


\2 


-P 2 kL- 


Note that Eqs. (2.4) and (2.5) differ from the analogous expressions in Ref. 1 in that, now, terms 
mSX appear due to swirl. The analysis later in Appendix A will show why this is the case. 
Additionally, the notation has been modified so that k replaces m as a subscript, where m, as 
before, is given by m = sB - kV. 


Standard vordcity waves will be represented here using the axial disturbance velocity 
component . All other information, such as the transverse component, will be derivable from 
this. This component is given by an expression analogous to that for pressure. 


$ = —-00 Jfc=— oort=l 


™AwsknU n (r)e L J . 


( 2 . 6 ) 


Here W = 3, and the normalization parameter cl, is the far-field speed of sound. The state 
vector, Awskn > represents the axial velocity mode amplitude u^ skn divided by m, and the U, f (r) are 
the duct radial mode functions associated with the vordcity waves. The normalization by m in 
A^ skn is needed to eliminate undesired division by 0 when m = 0 in some of the actuator disk 

formulas which are derived later. Note that the U„ (r) ’s are not the solutions of a particular 

differential equation associated with this analysis, but are an orthogonal basis set convenient for 
use here. 

Similar to the case for pressure waves, the axial velocity component in Eq. (2.6) is a linear 
combination of vorticity duct modes of the form 

U n (r)e Km ^^~ sBQl \ (2.7) 


where here the axial wave numbers for wave type W = 3 are given (Ref. 6) by 


y skn =■ 


m£l s — sBCl 
~U : 


( 2 . 8 ) 


which guarantees that these waves are convected with the mean flow. The vorticity wave radial 
duct modes are given by 


U n (r) = cos 


(n — l)7t (r - 
r D~ r H 


(2.9) 
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The modes in Eq. (2.7) are the same as those derived and discussed in Ref. 6. However, 
the notation has been changed to that used here and the form is slightly modified. Also, the 
modes have the additional factor U n (r) because they are three-dimensional rather than two- 
dimensional. The particular set of functions U n (r ) was chosen to have the properties expected 
physically for the vorticity modes; also because they provide a complete family of functions in the 
mathematical sense. This selection is similar to Goldstein’s choice in Ref. 8, p. 223. 

Before finally defining source vector coefficients and scattering coefficients, some 
observations about scattering should be made. As per the discussion in Ref. 3, upstream-going 
pressure waves will scatter into upstream- and downstream-going pressure waves and into 
downstream-going vorticity waves. Downstream-going pressure and vorticity waves will do the 
same. With regard to indices, the stator scatters input waves with indices k, n, into s, again, 
but not necessarily into h and n*. This scattering is viewed in the stationary coordinate system, 
which is our stator-fixed system. The opposite is true for the rotor in the stationary coordinate 
system; it scatters input waves with indices s„ h, n, into h again, but not necessarily into s, and h- 
Thus schematically we have 

Stator. Si, kj, tt j ) Sj, Jc, n, 

Rotor: s t , fc, n, -» s, k„ n. 

Note that mode scattering for the rotor can give both positive and negative values for the output 
harmonic s. 

In SOURCE3D, the source vector coefficients, B ^ skn , are the mode amplitudes produced 

by the rotor wake acting on the stator, i.e., they are the values of A^in Eqs. (2.1) or (2.6) for 
the output waves generated by the wake. The scattering coefficients are the mode amplitudes, 
S mr i skn s k i n i ’ produced by a unit modal input wave, i.e. one with A$. J;Jk . n . = 1. Alternatively, 

these coefficients are the ratios A^ skn /A^. Sjkj for the case where the input state vector 

coefficients are not necessarily unity but are arbitrary. Schematically, the scattering process can 
be represented by the notation 



Clearly, the indices must conform to the input/output requirements described above. 

The elements A ^ skn , B & shl , and can be used to define several vectors and 

matrices that are useful. The vector Ayy is made up of elements Ayr s ^ , where the indices ( s,k,n ) 
range over all combinations that the program has selected as active. is defined analogously. 
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pp pp. 

and is the collection of elements over all active combinations of indices 

( Si,kijii ) and (s,k,ri). As an illustration, consider the simple example where there are two 
harmonics, s = 1 and 2. Assume that the only active values of k and n for s = 1 are (k - 1, n = 1). 
For s = 2, assume the active values are (k = 1, n = 1) and ( k-2,n- 1). Take (P = 2, W = 1). 
Then we can write 


(a 2 1 

A 1,1,U 


*uu' 



* 

* 

a 2 

1 A l,2,l,l I 

II 


c 22 

S Sl2 ~ 

* 

22 

*^12,2,1,15,1,1 

c 22 

°12,2,U;2,24 

*— 1 
<N 


^2,2,1 
k J 


• 

c 22 

^12,2,2,15,1,1 

^12,2,2,155,1. 


where the dots in the matrix for S/f represent zero elements. Using the vectors and , 
we can define more precisely the vectors A p and B? used in Fig. 1.1. They are given by 





fa 

II 

a! 

Af 

% 

II 

A 

Bl 

*3 

< 


These give overall vectors 


A = < 

A 2 

B = < 

B 1 

B 2 ■ 


A 3 

k 4 


B 3 

K > 


An overall scattering matrix 5 is also defined, whose level of inclusion of elements corresponds to 
A and B. It contains dements for all the input and output modes of all the active waves at all the 
various interfaces. It is displayed in Refs. 3 and 5. Consult these references for further details on 
S as well as on other items in this section. 
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CHAPTER 3 


SOURCE VECTOR COEFFICIENTS AND SCATTERING COEFFICIENTS 
FOR THE STATOR AND ROTOR ALONE 


In this chapter, we obtain source vector coefficients and scattering coefficients for what 
will be called the “stator alone” and “rotor alone” cases. These are the cases where coefficients 
are derived for the stator or rotor before adding mode scattering by actuator disks. Scattering 
coefficients for the actuator disks will be presented in Section 4.1. Then, Section 4.2 will discuss 
a procedure for combining these quantities with the present source vector coefficients and 
scattering coefficients to produce the final values for the stator and rotor elements. 

The approach for finding source vector coefficients and scattering coefficients is the same 
as that in Meyer and Envia (Ref. 1) for obtaining acoustic wave pressure coefficients and as in 
Smith (Ref. 6) for determining vorticity waves. We must first derive the upwash on the stator or 
rotor coming from the rotor wake for source vectors or from unit modal amplitude standard 
waves for scattering coefficients. The upwash is the component of velocity of these inputs that is 
normal to the surface of the vanes or blades. Then the upwash is used to obtain the loading on 
the vanes or blades by solving Kernel function integral equations, hi Refs. 1 and 6, the stator and 
rotor are modeled as linear cascades of flat plates which are divided chordwise into “strips” for 
the calculation. The loads take the form of dipole distributions. These loads are then coupled to 
the acoustic or vorticity modes to give output modal amplitudes. 

This coupling is calculated using formulas for output modal amplitudes and output 
vorticity wave modal amplitudes that are derived in Appendices A and B. The formulas for 
pressure there differ somewhat from those in Ref. 1 because now there is swirl present between 
the rotor and stator. Also, the formulas for the vorticity wave case differ from those in Ref. 6 in 
that now we must modify the procedure to cover a three-dimensional rather than a two- 
dimensional situation. 


3.1 Upwash and Loading 

For the stator, the upwash w is evaluated in the stator-fixed system and, as discussed in 
Ref. 1, has the general form 


w _ _ £ w ^i(k s z-2^sB!V-sBOj) (3.1) 

5=— oe 

on vane v, where w, is the Fourier coefficient of w, k, is the chordwise gust wavenumber, z is the 
chordwise coordinate along the vane with origin at center-chord and extending from -b at the 
leading edge to +b at the trailing edge, and v is the vane index, v =0, - • • , V - 1. 
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To determine the upwash in the form required by Eq. (2.1), we must first determine the 
velocity of its source - either the rotor wake or standard input waves - in the direction of the 
upwash, i.e. normal to the stator. Then we must transform this result to chordwise coordinates. 
Details for obtaining w, and k s for the rotor wake are given in Ref. 1 and will not be repeated here. 
From Ref. 1, we have for the rotor wake that 


= —W Sin(a £ + d c ^ e ~^ sS ^ yRD+ysD ~ XsPAC t2DCL CL +x SPAC,E tana C£.H) / r ) 
i[^Z?(sina5+oosa5tana c ^)/r] 

X € 


(3.2) 


and 


k s = — (sina$ +cosa$tanac£,)- 

r 


In the above, W s is the Fourier coefficient of the wake velocity, W, downstream of the rotor, a cl 
is the rotor blade relative velocity flow angle, and xspac is the blade trailing edge to vane leading 
edge axial distance. The parameters as, yj®, ysD, and b have been defined previously; W„ a s , 
<Xcl, ypD, ysD, xspac, and b vary with radius r. These parameters are shown on the schematic in 
Fig. 2.3. There is a factor e ~ isBxsPAC - E tsma cL,H ,r ^duded in Eq. (3.2). This factor sets the phase 
in the circumferential direction so that the exponential term in Eq. (3.2) equals 1 at the stator 
leading edge at the hub. The subscript H in xspacs and q lcls denotes that their values are taken 
at the intersection of the stator leading edge with the hub. 

Note that SOURCE3D evaluates the wake Fourier coefficients semi-empirically using 
experimental data fitted to a Gaussian wake velocity profile. It allows several choices for these 
profiles and also allows use of models for hub and tip vortex flows. Further discussion of the 
wake and hub and tip vortices can be found in Refs. 1 , 9, and 10. 

To determine the upwash, when loading is from input waves rather than from the wake, 
we first obtain the upwash velocity components for each s. We start with either the standard 
pressure representation, Eq. (2.1), for acoustic waves or the standard axial velocity 
representation, Eq. (2.6), for vorticity waves. Using the momentum equations as in Appendix C, 
we can derive the axial, u, and transverse, v, components of acoustic velocity from Eq. (2.1). 
Similarly, using the mass conservation equation, again as in Appendix C, we can derive the 
transverse component for vorticity waves from the axial components appearing in Eq. (2.6). 
Once u and v are known, the upwash is found by resolving these components into a vector normal 
to the stator, using the relation 

w = u sin as - v cos cts. (3.4) 

The result is 
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™i 


w s =-^ 


Ps 


cosa 5 +y skjTii sina 5 
sBa + U s y skini -m£l s 


, , , sk m (*l - hsBQt] 

I 4*^ 


(3.5) 


for acoustic type waves, and 

c r • 1 TJ / \ ( X 1~ X - ,ry 

w s = -a„[ry ski7li cosa 5 - m, sma 5 Jf7„. (r)e ^ n . (3.6) 


for vorticity waves, where w s represents the 5 th harmonic of the upwash. 

For the standard waves used in deriving the above relation, we have taken x — x in 

s 

Eqs. (2.1) and (2.6) so that the input waves, at this stage, originate at an interface located at x , 
the axial position where the stator leading edge meets the hub. For convenience, the W subscripts 
have been omitted from the expressions, but it should be clear where they would have been 
added. The y sk . n . ’s are those appropriate to the particular wave types; m, is the input wave 

circumferential wavenumber. The variables p 5 and U s are, respectively, the mean density and the 
uniform axial velocity in the stator element region; the other parameters have appeared 
previously. 

We next change from stator coordinates, x, to the chordwise coordinate, z, using the 
transformation 


jcj =x s -x SD +zcosa s +bcosa s , 


(3.7) 


<t> = - 


ysp 

r 


z . 2tcv 5 fcana s 

h — smote 1 

r S V r 


(3.8) 


referring to the geometry of Fig. 2.3. Then we have 


w* = 


.Ezl 

Ps 


-y-cosa 5 +y skn sinas 




sBQ + U s y sk . n . -m,a 

‘Mr*,-!., x SD~ m <ysD /r Hmi /r)fc sin a 5 -Y .fc COS a s ] 


X e 


Ay^-n, 


(3.9) 


j'[(-^-sin a 5 ~7 skjn . cos a $ )z— 2 iw>sBIV—sBQj] 


X e 


for acoustic type waves, and 
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(3.10) 


w* = -a^ry^ cosa 5 -/n I sina s ][/ n .(r) 

JH Am X SD ~ m iysD ,r+ ( m i /r)bsina s -y A .„.bcosa s ) 


x e 


A S^ 


i[(— sina 5 -y sk . r! . cos a s )z-2tuvsB/V-sBQi] 
X e r 


for vorticity waves. 

By inspection of Eqs. (3.9) and (3.10), we now can write w, and k s . Because loading is 
from input waves rather than from the wake, then w s and k t depend on k and m in addition to s. 
To emphasize this dependence, we replace w, and k s with w s kjrl . and k 5 fc;n . . For input waves of 

acoustic type (W= 1, 2), we have 


w 




Ps 


"*»• 


cosa 5 +y skjn . ana 5 


sBQ + U s y sk . n . —mfii 

' 1-7 Am x sd ~ m iyso /r+ i m i tr)bsma s -y Afn .bcosa s ] 


Wmn<r) 


x e 


AskiTli 


— w s,k i n^sk i r i,- 1 


and for vorticity waves (W = 3), 

cosa s - to, sma$] 
= w s,k l n I Ajfc,n, ■ 


(3.11) 


(3.12) 


For all three wave types. 


= -Y costx s • 


(3.13) 


For the rotor, the upwash is evaluated in the rotor-fixed system where it has the general 


form 


v i(k k z~2tvvkV I B-kVClt) 

w = ~L w k,s ini e 


(3.14) 


for vane v, v = 0, • • • , B - 1. This is analogous to the representation in Eq. (3.1), only now sBQ. 
is replaced by kVQ., w k >J;n . and k ks . n . are the rotor counterparts of w s k ^ and k s k;Jl . , and z is 
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the chordwise coordinate along the blade with origin at the center-chord and extending from -bn 
at the l eading edge to +b R at the trailing edge. Because there are no wakes loading the rotor 
blades, there are only upwashes produced by standard wave type inputs. Otherwise, there would 
have been parameters w k and h also. 

We obtain the upwash in the form of Eq. (3.14) and the expressions for w k S;n , and k k s _ n . 

by employing a procedure parallel to that for the stator. We start with the same expressions as 
derived before for the standard wave components u and v in the stator-fixed system, then shift 
these expressions to the rotor-fixed system using the transformation <J> = + Qt , where 

temporarily <j> is the polar coordinate in the stator-fixed system. The exponential 
in the expression for the standard wave in the stator-fixed system then becomes 
_ e i[m t ^-(s i B-kV)Qj-s i BQj] _ ^{m^-kVD:) ^ rotor-fixed system, giving the 

factor e _ifcvn * appearing in Eq. (3.14). Using the resulting expressions for u and v, and 
proceeding as in the stator case from this point forward, we then obtain relations comparable to 
Eqs. (3.1 1)-(3.13). In proceeding in this fashion, however, we use the relation 

w = usin(X R -vcosa^ (3.15) 


instead of Eq. (3.4) and the coordinate transformation 

= x R - xjusd +zcosa* + b R cosa* 


, yRLED 2icvk b R sina^ 

(j) 1 — smoc % r 

r r V r 


(3.16) 

(3.17) 


in place of Eqs. (3.7) and (3.8), again referring to the geometry of Fig. 2.3. In Eq. (3.16), / is 
the axial location where the leading edge of the rotor meets the hub. 


The end result is that for input waves of acoustic type (W = 1, 2), we have 


m. 


w k,s ini = ■ 


.Ez 

Pi? I 


-cosa/f+y^ sina* 


N'm.n.W 


sBQ+U r y Sjbt . - m&. 

if— Y Sibti X RLED -f'hyRLED / r+(mj / r)b R sin a s -^ Sihl b s cosa R ] 

ft * * * 


x e 


V, (3-18) 


— w k,s l ru/‘ i s I foi l 


and, for vorticity waves (W = 3), 
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(3.19) 


w k,s i n i = -a~[n sfr cosa R - mi since* Jt/,. (r) 

v ^‘t-Y i; tei XRLED ~ m iyRLED ! r+ (" t f *r)b R since * -y *. to . fc* cosa 8 ] 

x e "Sjkni 

= ^ k.Sjn^Sjbij • 

Further, 


***», = y since* -y sM cosoe* 


(3.20) 


for all three wave types. The parameters a*, xrled, yw£D, and b^, used above, have been defined 
previously; they are functions of radius. The quantities p r and Ur are, respectively, the mean 
density and the uniform axial velocity in the rotor element region. Analogous to the situation for 
the stator, the input waves for the upwash calculation originate at an interface located at x R , the 
axial location where the rotor leading edge intersects the hub. 

Returning now to the stator, once the upwash is known, we can use the theory in Ref. 1 to 
find the loading. In particular, we determine /*(r,z), the s b harmonic of the elemental vane 
chordwise loading function, which is defined in Ref. 1; refer there for more detail. We mention 
here only that/,(r,z) is related to the s* harmonic of unsteady loading, A p s , on vane v = 0 through 
the relation = ApJpsU r sw„ where U r s is the fluid velocity relative to the stator vanes. 

After f s (r^) is determined, it can be used in formulas derived in Appendices A and B to obtain 
pressure and vorticity axial velocity modal amplitudes, which, in turn, give source vector 
coefficients and scattering coefficients. 

As per the discussion in Chapter 3 of Ref. 1, f,(r, z ) is determined as the solution to the 
integral equation 


«** — f K c (z-y)f s (r,y)A' (3.21) 

4 b 

where K c (z - y) is the cascade kernel function for the stator. For standard wave input rather than 
rotor wake input, k s and/* should be replaced by k skjn . and f Syk{n . - & is evaluated by using the 

following parameters (see Ref. 1): 

Pr : 


Reduced frequency: K = -sBQb / (P *c s ). 

Inter- vane gap: h = 2izr / V , 

Vane stagger angle: a s , 
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r: -~-D + M T M rS (rlr D )ana s /tf]. 

In the above, Cs is the no minal speed of sound in the stator element region, Mr is the rotor blade 
tip rotational Mach number, and M r s is the Mach number of stator relative flow, U^/cs. The 
negative sign is used in Eq. (3.21) because we need the loading that induces a normal velocity 
equal, but opposite in sign, to the wake or standard wave input 

For the rotor, the loading function for standard waves is rat ^ er than 

f s k (r,z) . If we proceed through the analogous steps for the rotor as previously described in 
Ref. 1 for the stator, we obtain the integral equation 

*‘ W = - I K c (z - y)f lm (r,y)%- , (3.22) 

b * 

that must be solved to determine . From the derivation of the Kernel in this case, it is 

easily seen that to evaluate K c for the rotor, the parameters above for the stator now must be 
changed to: 


p r : y/l -M? s , 


2 

Reduced frequency: K = -kVQ. b R / (P r c R ) , 

Inter-blade gap: h = 2jtr / B , 

Blade stagger angle: - a R , 

r : +?^[l-M T M rR (r/r D )sLna R /$ 2 r ]. 

B 

In the above, cr is the no minal speed of sound in the rotor element region and M t r is the 
Mach num ber of rotor relative flow, U t r Icr. The differences in T between the stator and rotor 
are due to the fact that for the stator the inter-vane phase angle is given by a s = -2 nsB / V , 
while for the rotor, its analog, the inter-blade phase angle, has the form g R = +2nkV / B (see 
Ref. 3 for discussion of inter-blade phase angles). 
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32 Source Vector Coefficients and Scattering Coefficients 


To find source vector coefficients and scattering coefficients from the loading, we use the 
formulas for acoustic output wave pressure modal amplitudes and vorticity output wave u/m 
modal amplitudes. The formulas and their derivations are given in Appendices A and B. These 
results are obtained by adapting the approaches found in Refs. 1 and 6, respectively, to the 
situation here. 

The source vector coefficients Bwsh, for acoustic output waves (W = 1, 2), as per the 
definition given in Section 2.2, are just the pressure modal amplitudes p Ws h, given by Eq. (A.24) 
when the input loading is from the rotor wake upwash. Similarly, the source vector coefficients 
B w ,kn for vorticity output waves (W = 3) are just the u/m modal amplitudes given by Eq. (B.12). 
Thus we have 


&Wskn ~ 


r H 

L 


r D b 

J Cwito( r ) w j( r )J fs( r 2)Dw s k n (r,z)dzdr for W = 1,2, 

-b 


(3.23) 


X fsin^Dwsbiin^dzWi for W = 3, 

l=i -b 


where the subscript W, which was implicit in the last section, is now shown explicitly. In the 
above, w f (r) is given by Eq. (3.2); Cwa*(r) and Dwsh>(r, z) by Eqs. (A.25) and (A.26) for W = 1, 2; 
and by Eqs. (B.13) and (B.14) for W = 3. The elemental vane chordwise loading function f s (r, z ) 
is deter mine d by Eq. (3.21). As discussed in Appendix B and illustrated in Fig. B.l, for the 
vorticity wave case, the annular duct is subdivided into L sub-annuli with widths A/j = (n - n.j), l 
= 1, ■ • • , L. At this stage, the interface for output waves is at X s , the axial location where the 
stator leading edge intersects the hub. 


To obtain scattering coefficients for the stator, we apply Eqs. (A.24) and (B.12) as above, 
only now the lo adin g is that provided by unit modal amplitude standard input waves with w sWj kin . 

given either by Eq. (3.11) or (3.12) with A w . skjn . = 1 . We have used subscripts i to denote input 

modal indices; since the input and output harmonics must be the same, no i has been used with s. 
We can write 


^WW^kn-.sk^n, 


r D b 

] C WsA7I (r)w 5W _ ^(r)J f sWiAril ( r ,z)D Wskn (r,z)dzdr for W = 1,2, 


X for W = 3, 

[l = i 


(3.24) 
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where W, = 1, 2, 3. Also, w sWjikjn .(r) is defined by Eq. (3.11) or (3.12), and CwihSj) and 
Dwdo>(r, z ) are the same as before. 

Thus far, the source vector coefficients given above are for output interfaces located 
where the leading edge of the stator meets the hub, i.e. at x s . The same is true regarding the input 
and output interfaces for scattering coefficients. Refer now to Fig. 3.1, which shows prior and 
new interface locations. For the stator alone case, we shift the interfaces so that downstream- 
going input waves originate from an interface at x 2 and upstream-going input waves from an 
interface located at the farthest downstream stator trailing edge coordinate with respect to radius, 
x e . Downstream-going output waves originate at interface x and upstream-going output waves at 
x . More generally, we designate the new input and output interfaces, respectively, as x ' and 


P P 

Because of the factor e -iyWjfa(Xl-JC } in Eqs. (2.1) and (2.6) for standard waves, it is 
easily seen that the modal amplitudes of output waves associated with source vector coefficients 
are changed by an amount 


E 


P 

Wskn 




(3.25) 


where Jwskn is the axial wave number in the stator element region. Scattering coefficients are the 
ratios of the output to the input of state vector coefficients, i.e. A Wskn / ; therefore, 

using reasoning similar to that above, scattering coefficients must be multiplied by the ratio of the 
output to input shift exponentials, which is the ratio 


fiPPi 

r WW i sk7i;sk i n i 


0 Wskn 


(x p -x s ) 


-n Widen ( xPi - xS ) 


(3.26) 


Thus Eqs. (3.23) and (3.24) become 


Bwskn - 


r D b 

Ewsbi j Cwskn( r ) w s ( r ) J fs( r > z ) ^Wskn ( r > ^ )*&^ r f or W = 1,2, 
r B -b 

L r 

EwsknY, forW = 3, 

1 = 1 -b 


(3.27) 


and 


- 21 - 



X 



-22- 


FIGURE 3.1 AXIAL INTERFACE COORDINATES ATTHE HUB 


cPPi 

^WW^kn^sk^ 


td b 

Fww^kn^n, J C Wsfet (r)iv^. M (r)J fsw^kjni ( r -,z)Dw s i cn (r,z)dzdr for W = 1,2, 

t b ~ b 

L b 

{CwrfnC//)^,^.^)/ f sW . <kjn .(r h z)D Wskn (r l ,z)dz}Ar l forW = 3, 
1 = 1 -b 


(3.28) 


where W, = 1, 2, 3. For cut-on waves, the F s of Eq. (3.26) merely amount to a phase shift 
However, for input or output waves that are cut off, the Fs include real exponentials. These 
quantities serve the essential role of keeping reflection and transmission coefficients small for 
decaying waves and avoiding large off-diagonal elements in the scattering matrix. 


Scattering coefficients for the rotor are derived analogously. Using Eqs. (A.45) and 
(B.15) of Appendices A and B, respectively, we have 

r D b R 

F$ iS kn;sM J C Wsbt(r)™kW h s ini (r) J fkW i ,s ! n l (r,z)D W skn( r ’Z)dzdr for W = 1,2, 
r H ~ b R 

L b R 

\FZ : skn;sMX ^Wsknin^kW^nSn) J fkW i ,s i n i ( r l,z)D Wskn (r l ,z)dz}Ar l 
2=1 -b R 

for W = 3, 


(3.29) 


c PPi -< 

'^WW i skn,s i kn i ~ < 


where W; = 1, 2, 3. Additionally, w kWj S . n Xr) is given by Eq. (3.18) or (3.19), and C Wskn (r) 
and Dwskn(ri,z) by Eqs. (A.46) and (A.47) or Eqs. (B.16) and (B.17). 


The term F^. skn . s kn . in Eq. (3.29) shifts the location of the input and output interfaces 
in analogous fashion to F^ skn . sk , n . for the stator case. It is given by 


pPPi 

"wWjSbViSjbii 


e ~ b (Wsh,( xP - xR ) 

-o 'Wi.jtaj (**-**) 

e 1,1 


(3.30) 


Without this factor, the interface for both input and output waves would be located where the 
rotor leading edge meets the hub, i.e. at x R . Instead, referring to Fig. 3.1, the output interface is at 
x 2 for downstream-going waves and at x\ the farthest upstream trailing edge location with 
respect to radius, for upstream-going waves. Input interfaces for upstream- and downstream- 
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going waves are the reverse. The variables x P ' and x p represent the locations x 2 and x‘ , 
depending on the context, and Jwskn is the axial wave number in the rotor element region. 

There is an additional item that must be taken care of before the standalone scattering 
coefficients are completely specified. For both the stator and rotor, we must add a term 

1 x Fww s kn y skn t0 transmission coefficients for those waves that are scattered into themselves, 

i.e. waves for which W t = W, s, = s, ki = k, tii = n. This term is added to account for 
transmission of the original wave through the stator or rotor. Though not appearing explicitly in 

Eqs. (3.28) and (3.29), this term will now be understood to be included where appropriate. The 
PP’ 

function Fwwskn;skn is defined either by Eq. (3.26) or (3.30), except that the subscript i is 
removed from all variables except P, in those definitions. 

Before closing, it should be noted that SOURCE3D employs non-dimensional geometry 
and performance parameters rather than the dimensional ones used above. This means that for the 
stator it utilizes lengths r D and bj and speeds cs, Us, and Urs, defined with respect to the 
stator element region. Here, br is the vane semi-chord at the tip. Using these parameters, the 
formulas defining the stator source vector coefficients and scattering coefficients can be written in 

terms of the dimensionless variables M s - Us/cs; Ocs = 2bi/r D ; c r = r^r^, y skn = y shl r D \ k skn 
= k skn r D ; Mj, the rotor rotational Mach number at the blade tip; M„ the swirl rotational Mach 
number at the blade tip; = K^rp; x = r/r D ; x s = x S i)l2bt, ys = ystJ^bf, x S or = x S p A c/2bf, 
and b = b/2br- Non-dimensionalization is carried through in the rotor case using analogous 
rotor parameters and dimensionless variables. 
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CHAPTER 4 


ADDITIONAL THEORY 


Two additional items of theory will be needed to complete the calculation of source vector 
coefficients and scattering coefficients: a discussion of flow turning using an actuator disk, and 
then a discussion of a combined element to couple either the stator and actuator disk or the rotor 
and actuator disk into one element. These topics will be covered in the sections below. 


4.1 Flow Turning at the Stator and Rotor 

The theory developed thus far provides swirl between the rotor and stator, but does not 
include the flow turning at the rotor and stator that must accompany it. In this section, results 
will be presented that extend the two-dimensional approach of Ref. 3 to the three-dimensional 
situation shown in Fig. 4.1. A first disk will turn the flow at the rotor leading edge and a second 
one will straighten it back again at the stator trailing edge. The jumps in mean flow properties are 
given by input to the prediction program. Then the associated jumps in unsteady properties, such 
as pressure, axial velocity, and tangential velocity, are determined using linearized versions of the 
equations of conservation of mass, transverse momentum, and axial momentum. 

With reference to Fig. 3.1, the actuator disks are located at the inlet interface x, which is 
the furthest upstream axial location with respect to radius of the rotor leading edge, and at the 
exit interface x e , which is the furthest downstream axial location with respect to radius of the 
stator t railin g edge. Input and output waves for the rotor actuator disk have their origins at the 
rotor upstream interface x 1 and at x'. Input and output waves for the stator actuator disk 
originate at x and at the stator downstream interface x . Note that these locations are not the 
final ones for rotor and stator input and output waves. These locations will be set by the 
combined rotor and stator elements discussed in the next section. 

Scattering coefficients at the actuator disks are based on the reflection and transmission of 
standard waves impinging on the disks. The harmonic indices s and k — and therefore m, stay the 
same for both input and output waves, so that actuator disk scattering is only among wave types 
(pressure and vorticity) and among radial mode orders. Similar to the two-dimensional situation 
(Ref. 3), there are three possible types of scattering. These interactions are shown in Fig. 4.2, 
where, for simplicity, only the subscript W, for W = 1, 2, 3, has been included with the state 
vector coefficient A; the indices P, s, k, and n, or n are implied but not shown. The parameters Ai 
and A 2 are pressure amplitudes, shown as solid arrows; the quantity A 3 represents the amplitude 
ulm for vorticity waves, shown as dashed arrows. In Fig. 4.2, incident pressure wave Aj gives 
transmitted pressure wave Ai, reflected pressure wave A 2 , and reflected vorticity wave A 3 . Inputs 
A 2 and Ai give similar results. Nine scattering coefficients will result from these interactions for 
each set of indices s, k, n h and n. The others will be zero. 
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FIGURE 4.1 FLOW TURNING AT ROTOR AND STATOR 
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THE THREE POSSIBLE TYPES OF SCATTERING INTERACTION 
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DASHED ARROWS VORTICITY WAVES. Aw, W=1 ,2,3, 

DENOTE STATE VECTOR AMPLITUDES.) 
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Scattering coefficients will first be derived for input and output wave origins at the 
actuator disk axial locations. Later, the wave origins will be shifted to the desired locations 
shown in Fig. 3.1. The details of the derivation of these results are given in Appendix C. The 
results in this Appendix hold for both the rotor and stator; only the input data for each will be 
different As discussed in Ref. 11, the actuator disk scattering coefficients are given by the 
elements of the matrix Ku>, derived there, and given by Eq. (C.58) here. This quantity satisfies the 
matrix equation 


A 0 ut — Ku, Ain, ( 4 . 1 ) 

where A^ and A out are vectors made up, respectively, of input and output state vector 
coefficients. They are given by 



\ 


Af 

II 

'a? 1 

\ a 2 I 

% A out =< 

A b 

a 2 r 


5* 


A3 


(4.2) 


where the elements A^ in the brackets above are the vectors defined at the end of Section 2.2. 
The quantity a denotes values on the upstream side of the actuator disk, i.e. those before the 
jump; b represents values on the downstream side, i.e. those after the jump. The elements of Ku, 

give the ratios A^ skn / A^ shlj . 

To obtain the final actuator disk scattering coefficients, the input/output wave interfaces 
must be shifted to x 1 and x' for the rotor disk, and to x 3 and x for the stator disk. This is done 
using the same approach as for the stator or rotor alone in Section 3.2, only now the just- 
mentioned axial locations are used. Thus each scattering coefficient given by Ku, must be 

DD 

multiplied by a factor given by 


for the rotor disk, and by 

******* 


(4.3) 


(4.4) 


P P. 

for the stator disk, where x and x ’ are the appropriate new input/output interface locations. 
Hence, the final actuator disk scattering coefficients are given by the relation 
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pp pp 

SwWjSbt-jkn; ~ F WW i slmjkn ; KioWW i skn,skn i » 


(4.5) 


where the witej are the individual elements of Ku,. 

As previously mentioned, before the coefficients above can be used in SOURCE3D, they 
must be non-dimensionalized. This is accomplished in the same way as discussed previously at 
the end of Section 3.2, using the same dimensionless variables as described there. The additional 

dimensionless variable wiAn = ^Wskn / c p » is also needed, where is defined in Appendix 

C and c P is the medium nominal speed of sound in the appropriate region. 


4.2 Combined Elements 

Thus far we have derived source vector coefficients and scattering coefficients for 
separate stator alone and actuator disk elements. We have also obtained scattering coefficients 
for separate rotor alone and actuator disk elements. Let us refer to the coordinates in Fig. 3.1. In 
the former case, the axial interfaces are at x and x e for the stator, and at x e and x 3 for the 
d i sk In the latter case, they are at and x 2 for the rotor, and at x 1 and x for the disk. 

However, SOURCE3D combines the stator and actuator disk into one element and does 
the same for the rotor and actuator disk. The combined stator/actuator disk element has axial 
interfaces at x 2 and x 3 ; the combined rotor/actuator disk element has interfaces at x 1 and x . 
Source vector coefficients and scattering coefficients for the combined elements are derived and 
presented in Appendix D. For the stator case, they are the elements of the source vector Bcomb 
and the scattering matrix Scorn given there, respectively, by Eqs. (D.19) and (D.18). For the 
rotor case, the scattering coefficients are the elements of the scattering matrix Scomb described 
there for the rotor. As should be obvious, the elements of Bcomb and Scomb for the stator 
depend on values of standalone stator coefficients and actuator disk coefficients. The same is true 
for Scomb for the rotor. 
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CHAPTERS 


CONCLUDING REMARKS 


This report has provided the analytical documentation for the SOURCE3D Rotor 
Wake/Stator Interaction Code. Source vector coefficients and scattering coefficients have been 
derived for the stator and rotor source elements, each in isolation, for use in the three-dimensional 
TFaNS (Theoretical Fan Noise Design/Prediction System). They have been obtained by 
extending the approach of the earlier V072 Rotor Wake/Stator Interaction Code. 

The V072 code uses only the stator as its noise source. The SOURCE3D program has 
added a rotor element, swirl flow and vorticity waves, actuator disks for flow turning, and 
combined stator/actuator disk and rotor/actuator disk elements. When used in TFaNS, the rotor 
element will allow frequency scattering and reflections from the rotor, neither of which V072 
provides. As mentioned previously, frequency scattering is the mechanism by which an original 
blade passing frequency (BPF) harmonic is scattered into other harmonics. The inclusion of the 
rotor element, swirl, and actuator disks accounts for mode trapping and amplification in TFaNS 
and explains features of experimental engine noise data not predicted by V072 (see Ref. 10). 
Mode trapping, as described earlier, is the phenomenon where the BPF mode grows between the 
stator and rotor blade rows through repeated reflections. This mode does not escape through the 
inlet or exit, but its energy is scattered into higher harmonics and is dispersed at two and three 
times BPF. 

It should be noted that the SOURCE3D program has been verified thoroughly vs. 
Hanson’s CUP2D two-dimensional coupled cascade noise prediction code using thin annulus test 
cases. Based on a modified advanced ducted propeller (ADP) geometry, the test model had 16 
blades, 22 vanes, and was run at a corrected rotor speed of 9,600 rpm. Hub radius te was 8.15 
inches and outer radius r D was 8.164 inches, giving an annulus that was sufficiently narrow for 
the case to be considered two-dimensional. For this situation, the two-dimensional predictions of 
CUP2D matched the three-dimensional predictions of SOURCE3D very closely. Comparisons 
were made for source vector coefficients and scattering coefficients for the stator and the rotor 
with and without swirl, with and without actuator disks, and for the actuator disks and the 
combined elements each in isolation. In all instances, the results were nearly identical. 

It should be pointed out that for two situations, SOURCE3D occasionally gives values of 
scattering coefficients which are unrealistically large. When this occurs, it is generally obvious to 
the user. The problems are caused by limitations in the original V072 model that carry over to 
SOURCE3D due either to two-dimensional load divergence or to duct resonance at three- 
dimensional cut on. Load divergence is a phenomenon that occurs in the two-dimensional 
unsteady aerodynamic theory when the normal distance between vanes or blades is half the free 
space acoustic wavelength at one or more radial locations. The problem at duct cut on arises 
from the factor k^, in the denominator of the three-dimensional Green’s function and in three- 
dimensional mode amplitude equations such as Eq. (A. 14). This factor, defined in Eq. (A. 10), 
goes to zero at three-dimensional cut on. In a fully consistent three-dimensional theory, this zero 
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lim it would be matched by a load distribution that also goes to zero at cut on. However, in the 
present model, our loading is based on two-dimensional strip aerodynamics. The two-dimensional 
loading does approach zero, but at the two-dimensional cut-on frequency. Because two- 
dimensional and three-dimensional cut-on frequencies are markedly different, computed pressure 
harmonics diverge at duct resonance. More will be said about these items in Ref. 12, which 
discusses three-dimensional scattering physics results using SOURCE3D. 

In conclusion, the SOURCE3D code adds the critical new ingredients of swirl, frequency 
scattering, and flow turning to the previous capability. This allows us to make predictions that are 
much more complete than before. 



APPENDIX A 


PRESSURE MODAL AMPLITUDES FOR A SYSTEM WITH SWIRL 


To obtain pressure modal amplitudes for the case where there is swirl between the rotor 
and stator, we apply an approach very s imilar to that in Meyer and Envia (Ref. 1) for the case 
where there is no swirl. There are separate formulas for the stator and rotor cases. As stated in 
Section 2.1, swirl is taken to be of solid body type ( V = £2,r) for the three-dimensional case here. 
It is implemented by letting the stator rotate at angular speed £2/ = -£2„ the negative of the swirl 
angular speed, and the rotor at angular speed Q (see Fig. 2.3). The difference, (£2/ - £ 22 ), is 
equal to £2, where £2 is the rotor angular speed. For both the stator and rotor, the analysis 
starts in the inertial system and ends, finally, in the stator-fixed coordinate system, which is where 
modal amplitudes must be evaluated. 


For the stator case, we begin with the Green’s function equation for unsteady pressure, as 
given by Eq. (41) of Ref. 13.* It is written as 


p(x,t) = -j J n(y) VG(x,y,t-T)Ap(y,x)dS(y)dx (A.1) 

— s 


in the inertial reference system. In this equation, x and y are, respectively, the field point and 
source point duct coordinates, x = (r,4>,X[) and y = (r',^',^) ; S is the airfoil surface area; n, 
the unit normal to this surface; Ap, the pressure loading on this surface; V, the gradient 
operator; and t and x, the field and source point times. The axial coordinate origin is taken to 
be at the location where the stator leading edge meets the hub, i.e. at x } = x s (see Fig. 3.1). For 
simplicity, the indices P and W have been omitted. It should be clear from context what they 
should be. G(x, y , t - x) is the Green’s function given by 


1 00 00 

G(x,y,t-x) = ]£ £ 




t ( r ) e »n(< 1>— 4>0 




g'Ty mn 


where 


Y mn(“) 


Ms 

fc s 


p 1 2 


sgn(y 1 -x 1 ) 


(A. 2) 
(A.3) 


and 


*mn(“) = 




/ \ 2 




(A.4) 


* This reference was the predecessor of Ref. 1 as documentation for V072. 
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Also, r = (tq - rtf )', to is radian frequency; cs is the speed of sound in the stator element 
region; Ms is the axial flow Mach number, Us/cs, where U s is the uniforrm axial velocity in the 
stator element region; and |3 = yjl-Ms ■ The other parameters have been defined previously. 

Let us shift the source point y to the stator-fixed frame using the transformation 

<j)'=<t>'+Q 2 T =<\>'-Q s x, (A. 5) 

where <j>' is the source point polar coordinate in the stator-fixed frame (see Fig. A.1). Also, we 
write Ap(y,x) in Fourier series as 


ApO>,t)= 2, A Ps(y) e sBQ * > ( a - 6 ) 

where y is takpn to be on the vane mid-surfaces S M , which have unit normals h(y) . Then, 
combining Eqs. (A.1), (A.2), (A.5), and (A.6), while rearranging, we can write p(x,t) as 


oo oo oo 


p{x,t)= 2 X X j 


1 


S— — °o m=— <»n=l 


4raTfc mn (co) 


'O' «*(©)* l~ ie * 


X | e*"-™*"*-* dxdS(y)da, 


(A-7) 




where Vj, is the gradient taken with respect to the y coordinates. 


Next we integrate over dx, then over 
where 8(x) is the Dirac delta function. Then 


d(a, while recalling that J e* 1 dx =2iz?>(x) , 
(0 = sBQ -m£l s and we have 


p(x,t) 


I II 


«> m=-»n= 1 


V m»(ry (m *" Wl ) e _i(sBa ‘ mQj ): 


1 

'Ifi'Krms 


X j Vmn( r, )”Cy) +7 ™ ,j) ' l) ]Ap i O)iiS(y), 

S M 


(A. 8) 


where 




FIGURE A.1 STATOR GEOMETRY FOR ACOUSTIC WAVE MODAL ANALYSIS 
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(A.9) 


Y mns =yrrm(sBQ-mn s ) 




sBQ.-mQ. s ^ 


M s 


±k mns 


< c s J 



and 


^nms=Km^BCl-n£l s ) 



(A. 10) 


Hie (+)-sign in Eq. (A.9) is associated with upstream-going waves (W= 1) and the (-)-sign with 
downstream-going waves (W = 2). 

Let us now shift the observer point y to the stator-fixed system where the modal 
coefficients are evaluated. We use the coordinate transformation 

<j) = 4> + = <{) — Q. s t , (A. 11) 


where <|) is the field point polar coordinate in the stator-fixed frame. Then Eq. (A.8) can be 
written as 




XIX 




i 

2ft'k mns 


X j Vmn( r ')w(j) V :) ,[e' ( '"♦■‘T— »)jAp a (y)dS(y) 

S M 


(A. 12) 


and p( x, t) has the form 

J>0M) = tH JWT-.C 

S=—oo 77i=— ooft=l 


In Eq. (A. 13), p ms is the modal amplitude coefficient we desire. It is given by 

Pmns — J ¥«n(0»(j) V y [e ,( "^' n ^ ) ]Ap,(y)ciS(y). 

^mns c 


(A. 13) 


(A. 14) 


This is the same as Eq. (4.29) of Ref. 1, except that now (sBQ - mil,) has replaced sBQ 
in the definitions of and The result of these transformations is that the axial 

wavenumber (and therefore cut-on frequency) given by Eqs. (A.9) and (A. 10) has been modified 
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to a form exactly analogous to that of the two-dimensional theory. The procedure retains the 
correct frequencies, sBQ . , in the stator frame. 

From this point on, we proceed exactly as in Ref. 1. Along with other relations, we use the 
expression 


n(y)-V, e l< ~ +Tmn,yi) j = i^cosa5+Y ffmj sina5^^ (A.15) 


the coordinate transformation 


yi = ~*sd +bcos<x s +z'coscl s , 

<t >' = (-ysD +bnna s +z'sina. s )/r' , 


(A. 16) 
(A. 17) 


and, the equality 


v-i 


A Pi = X Ps U rS w s( r )fs( r ’Z)e 
v=0 


-2 vkvsBIV 


(A. 18) 


where v is the vane index, -2nsB/V is the inter-blade phase angle for the stator, z is the 
source point chordwise coordinate along the vane, and the other parameters have been defined 
previously. For standard wave input, modal index subscripts are added to w,(r) and f s (r, z) in 
Chapter 3. 

In this process, we also evaluate a summation 

X e 2jtiv(m-,B)/v > (A. 19) 

v=0 

which equals V when 

m = sB - kV (A. 20) 


and, otherwise, equals zero. The result is Eq. (4.36) of Ref. 1, which gives /w as 

Pmns VwnWrS^W^cosa^+Y^sina^j 

r E 

b 

Xe i (-1mns x SD +m )>SD /r ) J ^( r z ) e 1 T(Ym« cosa S- m S ina 5 / '')(2 +b )]^ r 

-b 


(A.21) 
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Eq. (A.20) permits us to change the convention for the subscripts on /w, y™,, and k^,, 
and to write Eq. (A. 13) in the form 

Pw(x,,) = til PwstoV (A.22) 

k=-°° n=l s=— oo 

where jwskn and kwan are given by the right-hand sides of Eqs. (A.9) and (A. 10). Further, 
in Eq. (A.9) is now k W sh, because Ws have been made explicit Additionally, W = 1 or 2, 
because output waves in this section are pressure waves. With this new notation, Eq. (A.21) 
can be written, finally, as 


PWskn ~ 


PsV 


J V m „(r)17 rS w,(r)^cosa 5 +Y w , fat sina s j 

r H 


b 

Xe K-^w^ SD +my SD /r) J ^( r ^^'[(Y^cosa 5 -msma 5 /r)(z+fc)]^ r 

-b 


(A.23) 


Note that this modal pressure amplitude is the same as that in Ref. 1, except for one major 
difference. Because of swirl, (sBQ. - m£l s ) has replaced sBQ. in Jwskn and kwsim- Further, 
observe that w s and/, will depend also on input wave modal indices and n*, in addition to s, for 
the case of standard wave input, as discussed in Chapter 3. 

Eq. (A.23) can be written for later convenience as 

r D b 

Pwskn = J Cw,bi( r K( r )J f s (r,z)Dw sin (r,z)dzdr , (A.24) 

r H ~ b 

where 


<WO = „,<r(-cosa s +y sma , (A.25) 

D Whl (r,z) = e !l(YB ' jb,c ° sa ' s “ msinas/r ) (z+i ’) ] (A.26) 

For the rotor, we proceed similarly. Analogous to the stator case, the axial coordinate 
origin is taken at the point where the rotor leading edge meets the hub, i.e. at Xi = x R (see 
Fig. 3.1). We start with Eq. (A.1) and shift y = (r',<j>\;y 1 ) this time, however, from the inertial 
to the rotor-fixed frame. Here Eqs. (A. 5) and (A. 6) are replaced by 
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and 


ApOm)= X Ap* (y)e 


-kVQx 


(A. 27) 
(A.28) 


where <j)' is the source point polar coordinate in the rotor-fixed frame (see Fig. A.2). 
Using these relations, the analog of Eq. (A. 8) is 

p(x,t)= jr 2 I 


k=—<x> m=— oon=l 

X 


2/TC:, 


rrmk 


J VimOOwOO-V, e' ( m ^ +TmBkyi) Ap k (y)dS(y), 


(A.29) 




where 5 m now represents the rotor blade mid-surfaces, and and are given by Eqs. 
(A_3) and (A.4) with 


(0 = kVQ. + mill . 


(A.30) 


Similar to before, we shift the field point x from the inertial to the stator-fixed coordinate 
system, using that 

♦ -♦-a,*. (A.3i) 


where <|> is the field point polar coordinate in the stator-fixed system. Then the rotor counterpart 
of Eq. (A. 13) becomes 

= £it (A.32) 

fc=z — co TTl— — ooft=l 

where /w is the modal amphtude for the (fcVQ + mQj + harmonic. It is given by 

Pmr± I Vmn(On(y) V y p ( - w ^^>lAp t (y)d5(y). (A. 3 3) 

ZU s u L J 

We will see later that Eq. (A. 33) gives the modal amplitude we desire. 

Eq. (A.33) is analogous to Eq. (A.14). Thus we proceed exactly as for the stator case, 
only in place of Eqs. (A.15)-(A.18), we use 
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FIGURE A.2 ROTOR GEOMETRY FOR ACOUSTIC WAVE MODAL ANALYSIS 
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n(y) ■ V jy [e i(_wV+Y ™*^ ) ] = i^cosa* +y mnA sina R Je , ' ( ' mV+7 ^ l) , (A.34) 


>1 ~~ x rd +b R cos<x R +z'cosa R , 
ty' = (-y RD +b R sina R +z , sinai R )/ r' , 


B - 1 


A/>jfc = X PR U rR”k(r)fk(r,z)e 
v=0 


-UkvW/B 


(A.35) 

(A.36) 

(A.37) 


where z'is the source point chordwise coordinate along the blade. The exponent in Eq. (A.37) is 
different from that in Eq. (A. 18) because the inter-blade phase angle Or for the rotor is 
+2izkV/B.* This procedure gives 


_£j£_ 

PWm * 2n^ 


rD 

h 

r B 


mrt 


(r)U rR w k (r) I ycostx* +y ^^sina^J 


b R 

x e K- J 1v m *XRD +m tBD lr ) J Wmnk cosa fi -msina fi / r)(z+bg 

~ b R 


(A.38) 


where W’s have been added. 

In deriving Eq. (A38), we evaluate a s umm ation analogous to that in Eq. (A 19), i.e.. 


which equals B for 


I 


v=0 


2Ttiv(m+kV)IB 

& 5 


m = sB - kV 


(A39) 

(A40) 


and is zero, otherwise. Hence m must assume the same values as it did for the stator. 


Using Eq. (A40) and the fact that Q = Q/ + Cl s , we find that the time-dependent 
exponential in Eq. (A32) can be rewritten as 


g-i(kVQ.+mQ. l +m£}, s )t _ ^-i(fcV3Q+m£2)r 

_ e -i[kVQ+(sB-kV)Cl}t _ e ~isBQl 


(A.41) 


* The sign is the reverse of that for o s for the stator, because the relative motion between the 
rotor and stator is reversed. 
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Hence, we see that pwmnk is the modal amplitude we wish. Further, recalling Eq. (A.30), we have 
that 


G) = kVQ + m£2i = kV Q + m(Q. — Q, 

1 5 (A. 42) 

= kVQ + (sB - kV)Q. - n£l s = sBQ. - m£l s , 

so © is the same for the rotor as for the stator. Thus, referring to Eqs. (A.3) and (A.4), jmnk 
and kmnk are the same as for the stator, i.e. they are given by Eqs (A.9) and (A. 10), only Ms, cs, 
and P there must be replaced by their rotor counterparts. 

Changing the convention for the subscripts on pwmnk, Yhw, and as we did for the 
stator, we have, finally, that 


oo o© 


p w (x,t)= 2 X X PWsknV mn( r ) e,{m *^ WSh,Xl 


(A.43) 


k-~oo n=l s-- 


and 


_ P rB 

PWskn nT ^y 

2 TKbi 


J ^ nrn(r)U r RW k (r)y—cosa R +y Ws ^ana R j 

r H 

Xe i (-^Wsb, x RD +m yRD /r ) j yj t ( r?z ) e 4(T w *,oosa R -msina R /r)(z+fc i? )]^ r 


(A.44) 




In Chapter 3, subscripts s, and n, have been added to w* and j* to indicate that the upwash and 
associated loading are produced by input from standard waves that have these these modal 
indices. Similar to the case for the stator, pw S h> can be rewritten as 


r D b 

Pwskn = J Cwskn( r ) w k( r ) J fkiPZjDw^r^dzdr , 


(A.45) 


r E 


-b 


where 


Cwskn( r )- 


* . — ¥mn(q — cosa^ +Y Wskn 
2 Tk r ^ V r 


K-y Wskn x RD +m yRD /r ) 


(A.46) 


z) = e' [(YWtb,cosa/J_msinaR/r)(z+1 ’ fi)1 


(A. 47) 
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APPENDIX B 


VORTICITY WAVE MODAL AMPLITUDES 


For the present version of SOURCE3D, we use a simplified two-dimensional/three- 
dimensional approach to compute the axial velocity modal amplitudes for the vortidty waves. 
This is done by placing “very thin” sub-annuli at each radius r of the duct (Fig. B.la) and applying 
a two-dimensional Smith formula (Refs. 3 and 6) to the stator or rotor in each sub-annulus 
modeled as a two-dimensional cascade. This gives us the axial velocity at each r across the duct 
We match these values radially, on an (s, fc)-mode-by-mode basis, to values given by the three- 
dimensional expression for the vorticity wave axial velocity, which, at this point has unknown 
modal coefficients. By doing this, we obtain an equation involving a Fourier series of the radial 
vorticity modal functions U„(r) from Chapter 2. This equation can be inverted to give the 
unknown modal coefficients, which are the vorticity modal amplitudes we set out to dete rmine 
In the course of this process we discretize the equations by dividing the annular duct into L sub- 
annuli (Fig. B.lb); each sub-annulus has radius r t and thickness A r* = (r ; - r^j), / = 1, - • • , L. 

The Smith formula approximations are calculated in each of these sub-annuli. 

The two-dimensional axial velocity u(x,t) associated with a vorticity wave at a single 
frequency © for a two-dimensional cascade (Fig. B.2) can be deduced from S. N. Smith’s report 
(Ref. 6) to be 


u(x,t) = ] jr co S e + psine )Zo ^ o (B.l) 

0 S k=-— 

for the stator in the stator-fixed coordinate system. The notation is that of Ref. 3, where c is 
blade chord length; T , vortex strength representing loading; zo, vane chord coordinate; s, vane 
spacing; w 3 , a function given by Smith; a 3 , axial wavenumber, p, transverse wavenumber; co, 

angular frequency; and 0 is stagger angle. The reader is referred to Ref. 3 for further details. 
For a sub-annulus of fixed radius r, treated as a two-dimensional cascade, Eq. (B.l) can be 
rewritten in our notation with all the harmonics of blade passing frequency as 


■<(*,<)= X X 


(B.2) 





FIGURE B.2 TWO-DIMENSIONAL STATOR CASCADE 
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(B.3) 


The coordinate system is that shown in the “unrolled” geometry of Fig. B.3, with x specified by 
the cylindrical polar coordinates (r^jq) . The axes are attached to the stator at the leading edge 
at each radius = r location. As before, m = sB - kV. Also, J S kn = (mQs ~ sBQ.) / U $ , by 
virtue of Eq. (2.8); other variables have been defined previously. As mentioned in Appendix A, 
w, and f s will have additional subscripts, k and n„ for the case in Chapter 3, where the input is 
standard waves. 

The three-dimensional axial velocity, u(x, t), will have the same form as given in Chapter 2 
for standard vorticity type waves. We write it as 


■.(*.») = «- £ t t^W***™'*"^ (B.4) 

k=—a° s=—°° n=l 

in the stator-fixed coordinate system. Here, x = (r, <j> , x 2 ), where x 2 has its origin at the 
lifting edge of the stator at the hub (see Fig. B.3), which is where output waves originate for 
stator results in this section. The coefficients are the modal amplitudes we wish to 

determine. The inclusion of m in these coefficients is required to eliminate division by 0 when 
m = 0 in some of the formulas below. Subscripts W, which are equal to 3 for this case, have 
been omitted, but are implied. They will be added later. 

We rewrite Eq. (B.4) in terms of the coordinates for x , by applying the tranformation, 

x l =x 1 -x SD , (B.5) 


a T ysp 
<J) = <j> , 

r 


(B.6) 


to Eq. (B.4). These relations were obtained using Fig. B.3. Eq. (B.4) then gives 


u (x jt )= jr jr jjr f “^Ly (_my - sc/r+7 ^^ D) t/„ w ^ 


jb=— » ln=l 


-sBQi) 


(B.7) 
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FIGURE B.3 STATOR GEOMETRY FOR VORTICITY WAVE MODAL ANALYSIS 
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To determine the coefficients ( unJm ), we require the values of u specified by Eqs. (B.2) 
and (B.7) to match at each radius r. Because the exponential terms in these relations are the 
same, this means we must have 


“afc(r) = X 

n = 1 



e ‘(~ m ysD^ r+ y skn x SD^JJ (j-'j 


(B.8) 


or, alternatively, that 


u st (rje'' m,SClS ’~' Uk ' so ’ 




= X 

n=l 


u skrt 

m 


Icos^M 
) L r D~ r H . 


(B.9) 


because the exponential in Eq. (B.8) is independent of n. The right-hand side is a Fourier series 
with coefficients (u skn / m ) > so that, in essence, the two-dimensional r-dependent tranverse 

velocities have been resolved into radial modes. 


Eq. (B.9) can be inverted and rearranged to give 

= — - — f u sk (r)e i{mysD,r ^ *” XsD) U n (r)dr 
V m ) r D -r D £ 


(B.10) 


In practice, u sk (r) is evaluated only at the radii r u l - 1, • • • , L, shown in Fig. Bib. Hence Eq. 
(B.10) must be rewritten in discrete form. It becomes 

fes-l = — — Y u sk ( ri )e i{ - mysDlri ^ sbtXsD) U n {r l )Ar z , (B.l 1) 

\ m ) r D -r D /=1 

where Ar ; = - r / _ 1 . 


Finall y by recalling Eq. (B.3) and including Ws, the above can be written more 
conveniently in the form 



b 

^'Wskn( r O yv s( r l) | fs( r hZ)Dwskn( r bZ)dz 

-b 



(B.12) 


where 



sBQri 


C Wskn( r O~~ 


VU c 


u n {n) 


U, 


rS 


Amy SDln-^Wskn^SD) 


n( r D~ r H)n u rS 


'sBQnv 2 J smr A . 


u 


rS 


) +m ~VW” sina5 


(B.13) 


Dw skn (r h z) = g'[^w s far<»sa 5 -msina 5 /nXz+<')] 


(B.14) 


Because output waves in this section are vorticity waves, W = 3. 


The modal amplitudes for the rotor are dete rmine d by a similar process. We make 
changes in Eq. (B. 1) .s imil ar to those in going from stator to rotor in Ref. 3. We also change from 
w s (ri) and f s (ri, z) to w*(rO and fk(n, z), much as in Appendix A. In Chapter 3, subscripts s, 
and rii have been added to w* and /*, specifying the modal indices, because the input is from 
standard waves. Proceeding in parallel fashion to the stator case, we find that u(x, t) is still 
given by Eq. (B.4), only now xj has its origin where the leading edge of the rotor, rather than of 
the stator, meets the hub, so that output waves originate from this location. Further, Eqs. (B.12)- 
(B.14) are replaced by the relations 


Cwsknin)-- 


( u Wskn ''j _ 

f 

L 

=x < 

l m J 

/=1 

s 

bu r 

U n in) 

ft( r D ~ r Il) r l 

U rR 


Cwskn( r i) w k(n) J fkin^Uwshtin^dz 




kVQr, 


\& r h 


(B.15) 


u 


rR 


kVQrj 


U 


rR 




{ kVOrA 

-^n a * 


_ e ‘( m yRD lr l^r Wskn x RD 1 ^ (B.16) 


and 


Uwkn( r h z ) = e (7 WlJk " 00 S “ s _m Sinas /r ' X Z+i>R , 


(B.17) 


so that Eq. (B.15) now gives the rotor modal amplitudes we desire. 


- 48 - 



APPENDIX C 


SCATTERING COEFFICIENTS FOR ACTUATOR DISKS 


This section extends the two-dimensional actuator disk theory of Ref. 3 to three 
dimpmsinns The extension here was first developed in private notes by D. B. Hanson and then 
developed further by the present author. As discussed in Chapter 4, there are two actuator disks. 
The first is at the le adin g edge of the rotor where it turns the mean axial flow to add swirl 
between the rotor and stator. The second is at the trailing edge of the stator where it returns the 
flow to the axial direction (see Fig. 4.1). The theory here applies to either of the two disks. The 
configuration used is shown in Fig. C.l. For simplicity, the actuator disk is located axially at 
xj = 0, which is where input and output waves also have their origins. The procedure for shifting 
the input and output wave interfaces is discussed in Section 4.1. Region a in Fig. C.l will refer to 
the region on the upstream side of the actuator disk, region b to that on the downstream side. 

We will derive scattering coefficients for the three types of wave interactions shown in 
Fig. 4.2 and discussed in Section 4.1. To do this, we will use the jump conditions provided by 
satisfying linearized mass, axial momentum, and transverse momentum conservation requirements 
across the disks We do not include the effect of a radial component w. Radial equations are 
easily derived but are not part of the current formulation. The technique used is to write the flow 
variables as the sum of a mean part and a perturbation part (assumed to be small). Expressions 
for the conserved quantities on each side of a disk are equated. Mean flow quantities are 
determined and specified from a separate analysis. Jumps in the mean quantities at the disk 
necessitate jumps in the perturbation quantities according to the analysis of this section. The 
conservation requirements will be applied one (s, fc)-mode at a time; harmonic indices s and k 
(and therefore m) will be fixed across the disk, because scattering occurs only on radial indices n 
and wave type. First, though, we will need formulas for all the acoustic wave and vorticity wave 
unsteady velocity components in terms of state vector amplitudes. 

Velocity Components 

Acoustic pressure can be written in standard wave form (see Eq. 2.1) as 

pZto.r.M = p.f it , (C.l) 

s=-oo jfc=— « n=l 

where P will refer to either region a or b, and W will equal either 1 or 2. Superscript P is not 
needed here for ^^(r) , nor in the case later for U„(r), so it is omitted from these quantities. 

We seek expressions for the associated axial and tranverse velocities in terms of the 
parameters A ^ skn . We write these velocities in the form 
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ACTUATOR 




DISK 


DENSITY : 

Pa 



pb 

C h 

SOUND SPEED: 
AXIAL VELOCITY: 

C 3 

u| 

v a 

v a 




TRANSVERSE VELOCITY: 



v b b 


REGION a X 1= 0 REGION b X 1 


FIGURE C.1 CONFIGURATION FOR THE ACTUATOR DISK 
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= a. i it , (C.2) 

j=-oofc=-oon=l 

v£(*„r,<M> = a. t ft . (C-3) 

j=~oo fe=— oon=l 

and determine the unkn own coefficients, and v£r jfcn , through use of the momentum 

equation 


-+UVu = -—Vp, (C.4) 

dt p 0 


where u = (u, v, w) is the perturbation velocity; U = (t/, V, 0), the mean velocity; po, the 
n ominal fluid density; and p, the acoustic pressure. 


Eq. (C.4) can be written as 


du 

dt 


d_ 

dxi 


ZL +(t /^ +Q « )B = __LVp, 
«t> Po 


(C.5) 


3 15 

if we recall that = . The (-)-sign in the preceding relation is a result of the fact that 

dx 2 rdty 

x 2 and <[> are positive in opposite directions (Fig. 2.3). We have also used V = -Q,r. 
Eq. (C.5) gives the axial momentum equation 


— +t/— + q,^ — 
dt dx\ o<|> po 

and the transverse momentum equation 

dv TT dv _ dv 1 1 dp 

— +U — +£ 2 S — = -r~ 

d t dxi p 0 r 3<J> 


(C.6) 


(C.7) 


If we substitute Eqs. (C.l) and (C.2) into Eq. (C.6) on a mode-by-mode basis, while 
replacing U and p 0 with U P and p P , we find that 


p 

u Wskn 


P 

Poo Y Wsbi 
1 P y 

^ooP p A, Wjfet 


(C.8) 


where 


A-Wsfen - -sBQ, - U p y Wskn + m ^s • 


(C.9) 



Similarly, using Eqs. (C.l) and (C.3) in Eq. (C.7), we see that 


V P - P 

v Wskn ~ 


m 


P r ^Wskn 


(C.10) 


Thus we have 


s=-*ok=-oon=l P P'^Wsbi 


CC.11) 


and 


oo oo 


V w( X \’ r >ty>t) := Poe 111 


m 


PAW 


Aw l hW m ^.r)e K '*'’ ^ *• x '- ,Ba ‘ , . (C.12) 


To obtain comparable results for vorticity waves, we write the axial velocity in standard 
wave form [see Eq. (2.6)] as 


M = o_ 111 


(C.13) 


s =- oo k=-*x> n=l 


where W now is 3. 

We seek an expression for the transverse velocity in terms of the state vector coefficients 
Aw skn . To do this, we write the transverse velocity as 

vjta,r,*,0-«.2 X X . (C.14) 

j=-oo £=-oo n=l 

and find v^r skn through use of the continuity equation 

V« = 0, 

which must be satisfied by vorticity wave velocity components. This gives 

du _ 1 3v 

r 3q> ’ 

since the radial component for this wave type is zero. Substituting Eqs. (C.13) and (C.14) mode- 
by-mode into Eq. (C.16), we then have 


(C.15) 


(C.16) 
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(C.17) 


P _ P 
v Wskn ~ r Y Wskn • 


Thus we can rewrite Eq. (C.14) as 


v£(*,.r. ♦,*) = «.£ (C.18) 


s~-oo /:=— c>o n=l 


Conservation of Mass 

As shown in Ref. 14, the first-order equation for conservation of mass can be written as 


P a^a 


U a 

+ — Pa 

ci 


— ~Pa = P b u b 


2 Pb ' 

Cb 


(C.19) 


where a and b are the same as before; the mean flow parameters p a , p*, U a , U b , are considered 
known; and items with tildas represent disturbance, i.e. fluctuating components. Eq. (C.19) is 
applied at the interface between regions a and b, where xj = 0. 

Let us define 


C»=-^(PA+%A.>. ( C ' 20 > 

Poo C a 

c b =— (pA+%ft). < c - 21 > 

Poo C b 


where — is used for non-dimensionalization. Substituting values for p a and u a using Eqs. 

Poo 

(C.l), (C.ll), and (C.13), we have 


Q=^-[ Pa X ^^4** ^ (r) + X + 

boo Si P« X 1 sbt 


a ri= 1 

M„ 


PaX — + — P-X not ( r ) + P <P- X ’"A&wMOL 

„•=! Pa ^Iskn C a n'=l n'=l 

(C.22) 


for each fixed set of indices (s, k). The radial index n’ is used here to differentiate between this 
index and a second radial index n, which will be introduced later.* The factor e l ( m 4 , ~ sBa *'> , which 
is part of the disturbance components, is omitted from Eq. (C.22), but implied. 

* Note that, in the code, we have used the notation n; rather than ri. 
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We approximate C a , using finite sums, and write it more compactly as 


N\ _ N 2 _ N 3 __ 

Q = X rrm ( r ) A \sbi + X mri ( r )^2sbi + X ^3sht'^n( r )-^3sbi » (C.23) 

n=l n’=l n’=l 

where, after simplification, we have 


/-■fl ^ooY 1 jfcrt' i ^oa^g 


-Ukri 




l$An’ 


<£*»• = 


^ooT 2^' | gJgg 

c a 


ir* - 
^3skri 


aip a m 

P~ 



(C.24) 

(C.25) 

(C.26) 


In going from the first to the second line in Eq. (C.26), we have used the relationship 

p = 1.4p/c 2 


(C.27) 


for isentropic flow of air. There exist relations analogous to Eqs. (C.23)-(C.26) for Cb, C ^ skri , 
Clskn' ’ ^ C% shi . 


From a practical standpoint, we cannot satisfy Eq. (C.19) at all r, so instead we satisfy it in 

a Galerkin sense. We set C a - Cb, multiply by — f , and integrate over dr from r H 

Sin 


to r D , where 


r D 

Sin = | Vlm( r )rdr . 

r B 


(C.28) 


We do this for each n, n = 1, • • • , Nj t to obtain the set of equations. 
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«i n 2 n 3 

X - X - X 

n'=l n'=l n'=l 




W» 


AT, 


- X Clsknn' A Lkri ~ X ^2sknn^2sbi ~ X ^Sskmi A 3skfi > 
n=l n'=l n =1 

n = !,•••, 7^. 


(C.29) 


In writing this expression, we have placed the input A’ s, as shown in Fig. 4.2, on the right and the 
output ones on the left The parameter g ln in Eq. (C.28), and subsequent g’s, will normalize 
diagonal elements to 1 in matrices that will be defined later. In Eq. (C.29), 


riP c /~*P 

^ Isknri ~~°n,ri'-'lsbi > 

(C.30) 

^2 sknri ^ n y ri ^2 skri * 

(C.31) 

C?>shm ~~ ^\n y ri ^3sbi 

(C.32) 


for P = a orb. Further, 8 n „- is the Kronecker delta (1 when n = ri, 0 otherwise), and 


= — j U n - (r>|( (C.33) 

Eq. (C.29) will be used later when we set up an overall matrix equation. 

Conservation of Axial Momentum 

The first-order unsteady axial momentum equation at the actuator disk is given by 

(1 + M^)p a +2p a U a u a = (1 + M%)p b +2p b U b u b (C.34) 

(see Ref. 14). Thus we require that 

F a = F b , (C.35) 

where 

F a =—[(l + Ml)p a +2p a U a U a ], (C.36) 

Poo 

and Fb is the same, only with b’s in place of a’s. The factor 1/p., is used for normalization. 
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Substituting standard wave expressions into F a and Fb, truncating the sums as in 
Eq. (C.23), multiplying Eq. (C.35) by — — V|r mn ( r ) r for each n = 1, • • • , N 2 , and proceeding 

£in 

similarly to the conservation of mass case, we obtain the set of equations. 


Ni N 2 Ni 

X Kskmi A lOri ~ X F 2shm A 2sbi ~ X ^tsbm A 3sbi 


ri= 1 


n= 1 


n’=l 


AT, 


N, 


N* 


- X F lsknri A lskn ~ X F 2sbm' A 2sbi ~ X F 3sbm' A tshi » 

n=l rt'=l n'=l 

n = 1 ^ 2 . 


(C.37) 


This is the second set of equations needed for the final matrix equation. In Eq. (C.37), 


fu™ =6„.[(l + M^) + 2c P M f |^], 

^1 shi 

(C.38) 

*L« =5„,„.[(1 + M|)+2c p M p |^], 

k 2 shi 

(C.39) 

\PooJ\Cp J 

(C.40) 


forP = a or b. 

Conservation of Transverse Momentum 

The first-order unsteady transverse momentum equation is derived by retaining only the 
first-order terms from Eq. (D-23) of Ref. 14. It is given by 


P a U a^a + P a V a“a + P a U a V a ~ P b U b^b + P b V b“b + P b U b V b ■ (C.41) 


Since p = pc 2 , V a = Cl sa r, and V b = where and Cl,b are swirl angular velocities in 
regions a and b, this can be written as 


PaU a V a + Pa&s a ™ a +M a 


Q r „r 


Pa ~ P bUbh + P b&sb™b +M b 


-a J 




(C.42) 


Let us set 


G* 




fQ rl 


P a ^ a v a P a ^$<2 ^ a j 




v. c a J 


(C.43) 
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and define Gb analogously. Then Eq. (C.42) is satisfied when 

G a = G b . (C.44) 

We proceed as previously, substituting standard wave expressions in Eq. (C.44), truncating the 
sums, multiplying by ~^—U n {r) for each n, n = 1, • • • , N 3 , and integrating over dr from r B to 

82 n r 

r D , where 

« 2 , = 1 Ul(r)dr- < C - 45 > 

r B 

We obtain the result 



n'=l ri= 1 ri— 1 


n 1,* * 


(C.46) 


where 


p 

G\shm ~ -82 n,n mc P , P ? +5 3 n,n'KpCp , P > 

(C.47) 


^IjW "'lsht' 


P 

~P C- Mp R Hf r 7 2shi 

G 2sknri - 2*. n’ mc P ^ /> " + ° 3n t n‘ . p 

(C.48) 


A lsbi A'lsbt 


G 3sknri ^ n.,n’ 

+ 1.4— 

n /» — n 1 

(C.49) 


for P = a or b, and 


6 2n,n' = — J -TV m' «' (r)U n (r)dr , 

S2n i ' 

8 3 n,n=— Um'n'(r)U n (r)dr. 

* 


(C.50) 


(C.51) 


- 57 - 



In the above, M lP = Q. s prtJcp, the swirl Mach number at the outer duct radius. 

Matrix Equation 

Let Cj be the matrix consisting of all elements Cf ifetn . for all allowable s, k, n, and n,. 

t /ifc ^">3 fta fib t?cl r?b rfl 17 £ pfl pt fia fib f*a 

Lret Cj , 1 ^2 , tj , C3 , r; , ’ **2 * ^3 » ^3 > ^ , Cr^ > ^2 » ^3 » 

and G\ be defined analogously. Then Eqs. (C.29), (C.37), and (C.46) can be combined and 
written in matrix form as 


where 


K out A ou t K in A in > 


(C.52) 


K out 


Cl 

Ff 

G1 


-c\ 

-F2 

-G\ 


-c\ 

-n 

-G\ 


(C.53) 


Kin 


ci 


G\ 


-C°2 -Cl 


Fj -F? -Ff 


-G% 


-Gf 


(C.54) 


A i 

a|K 

A b 

A3 


(C.55) 


A» = 


A b 

A 1 

A| 

A| 


(C.56) 


The vectors A out and A u, consist, respectively, of output and input state vectors of type a£ • 
the vectors a£ have been defined earlier in Section 2.2. 


Eq. (C.52) can be solved to give 


Aon/ KgtgKfaAfa 

= K io^in’ 


(C.57) 


where 




(C.58) 


- 58 - 



Clearly the elements of K w give the ratios of the output to the input state vector coefficients. By 
the definition of scattering coefficients, this means that Ku> is the scattering matrix for the 
actuator disk This result is for input and output waves originating at the actuator disk. After 
these locations are shifted axially, in Section 4.1, we obtain the final matrix. 
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APPENDIX D 


COMBINING STATOR OR ROTOR WITH AN ACTUATOR DISK 


The source vector coefficients and scattering coefficients for the stator alone and 
scattering coefficients for the actuator disk have been derived separately for the configurations 

shown in Fig. D.la. In this figure, each represents the vector made up of all components 

p 

A-wskn °f type W at interface P over all allowable values of the indices s, k, n. The quantities 

are defined analogously for source vectors. See Section 2.2 for a more detailed description of 
these vectors. In this appendix, we derive the scattering coefficients and source vector 
coefficients for the combined stator/actuator disk element shown in Fig. D.lb. Afterwards we 
will discuss a combined element for the rotor. Note that the indices for the interface locations in 
Fig. D.l are kept general, i.e. a, b, and c, so that similar notation can be used later for the 
rotor/actuator disk situation. Then rotor results will be analogous to stator results. Referring to 
Fig. 3.1, we see that for the stator/actuator disk case, x a = x 2 , x b = x e , and = x 3 . The 
analysis here was developed jointly with D. Hanson. 

Referring to Fig. D.la, we can write for the stator 

A" = SM + S“ Aj + S?f A,‘ + B * , (D.l) 

A\ = S&A? -l-S^Af + S“a ; * +e|, (D.2) 

Aj=s£A2+S&Af+S;}A?+Bj, (D.3) 

where we have equated standard waves on the left with their component parts on the right 
Exponential and other factors, that would be present in the standard wave representations, are 

omitted because they cancel from both sides. In the above, the quantities 5^. , defined in 
Section 2.2, denote the scattering matrices made up of all the scattering coefficient elements 

pp m 

SwWjSknjkjn, for stator alone whose indices fall in the range of values allowed. The ’s 

on the left of the equations are for output waves; those on the right are for input ones. 

For the actuator disk, we have similarly 


A? = SgA} +SgAj +S{|AJ , 

(P-4) 

^2 = $ 22^2 +$23 A? +^22^2 » 

(D.5) 

A|=SgAl+SgA‘+S?,AJ. 

(P.6) 
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(B) STATOR AND ACTUATOR DISK AS ONE ELEMENT 


FIGURE D.1 COMBINING THE STATOR ALONE AND ACTUATOR DISK INTO ONE ELEMENT 
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pp pp 

where the scattering coefficient elements in the scattering matrices above are 

those for the actuator disk. Again, the A^ ’s on the left are for output waves, those on the right 
for input ones. 


Eqs. (D.1)-(D.6) can be written in block matrix form as 




‘ST 

s ab 

• 


<3 


B* 


> — 

s aa 

s bb 

s bc 


A b 

> +< 

B b 

A c 

out 


s cb 

s cc 


A c 

W J 

in 

B c 


where 


CD-7) 


Af 

4 

B P = 

B? 

b! 

S™ = 

r*aa c aa " 
* *12 *13 

, s ab = 

sf, ■ ■ 

? 

A/ 


B? 


. 


• . 






r abb 

*12 

i 




• 

c l 

II 

r*ba 
* *22 

0 >ba 

*23 

% 

II 

Q bb 

*21 

• 


II 

ii 


5 


c ba 

*32 

rtba 

*33 . 


,00 

• 



• 

* 


(D.8) 




- 


~ 

- 

s cb = 

aCb 

*22 

acb 

*23 

5“ = 

aCC 

*21 

. 


Qcb 
J *32 

r*cZ> 
*33 _ 


Sfi 

• • 


In the matrices above, the dots denote zero matrices. The notations { } fn and { } ou , in 
Eq. (D.7) designate, respectively, input and output waves. 

Eq. (D.7) can be rewritten in the alternate matrix form as 

AL=SrAl+S**i t +B°, (D.9) 

=S ba A? n + S bb A l + S* C A? +B b , (D.10) 

<«t = S cb AZ,+S ec A? t +B c , (D.ll) 

where we have changed the A p ’s from Eq. (D.7) to A? ’s and A^ M ’s to indicate whether 
they came, respectively, from the input or the output vectors there. 


- 62 - 



To obtain the scattering matrix and source vector for the combined element case of 
Fig. D. lb, we must eliminate the A? n ’ s and ’s from Eqs. (D.9)-(D.ll). To do this, note 

that, for the combined element, the output waves in Fig. D. la in the region between the stator and 
the actuator disk must be the same as the input waves going into the actuator disk there. 
Similarly, the output waves traveling from the actuator disk to the stator must be the same as the 

input waves entering the stator there. Hence, we require Af n = . Let us set 

A 6 =A} n = A b out , then Eq. (D.10) gives that 


(7 - S bh )A b = S ba A? n + S bc A? n + B b , 
so 

A b = ES ba A? n + ES bc A^ + EB b , 

where 

E = (7 - S bb ) _1 . 


(D.12) 

(D.13) 

(D.14) 


Using Eq. (D.13) in Eqs. (D.9) and (D.ll), then 

^ = (S“ + S^ESh* )A? + S^ES^Al + S^EBh +B a , (D. 1 5) 

= S cb ES ba A? n + (S cc + S cb ES bc )A? + S cb EB b +B C . (D. 1 6) 


Hence, 



[S M +s ab ES ba ] 
[S^ES*"] 


[ S ^ES bc ] 1 A a 1 \{B a + S^EE* } j 

[ S « +sC i £S * c ]J| A cj^ [{B c +5^ EB b }\ 


(DM) 


T akin g the two state vectors above as A out and Am, and setting 


> COMB ~ 


[S™ +S ab ES ba ] [S^ES^] 
[ S cb ES ba ] [S cc + S^ES^] 


(DM) 


_\{B a +S ab EB b )\ 

comb-^BC +s cb E B b }) 


we have, finally, 


Kui - S CO MB An + B COMB ■ 


(D.19) 


(D.20) 


When Bcomb = 0 , we see that the elements of Scomb provide the ratios of the output to the 
input state vector coefficients. Therefore, by the definition of scattering coefficients, Scomb is the 
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scattering matrix for the combined element Also, clearly, B C omb is the source vector for the 
combined element 


The derivation for the rotor case is similar. In place of Fig. D.l, we now use Fig. D.2 to 


obtain the equations 

Aj = S I 2 A 2 +S™A% +S$Aj , (D.21) 

A\ = S 22 A 2 + S^Af + S\\A\ , (D.22) 

A 3 = $32 A% + $33 A* + S 3 jAj , (D.23) 

A] = S b j b 2 A b 2 + S bb A b +S%A C j, (D.24) 

A 2 = S c 22 A 2 +Sf 3 A h 3 +Sf$A 2 c , (D.25) 

A 3 = S cb 2 A b 2 + S cb 3 A\ + SfjAj . (D.26) 


These equations have the same form as Eqs. (D.1)-(D.6) have for the stator case, only source 
vector terms are no longer present and x a = x 1 , x h = x, and x c = x 2 . This change in x a , x b , 
and x c means that rotor alone coefficients replace the actuator disk coefficients in Eqs. (D.l)- 
(D.6), and actuator disk coefficients replace the stator alone ones. 

Proceeding from this point on, exactly as before, we now obtain 


Aout - ScombAui (D.27) 

in place of Eq. (D.20), where Scorn is given by Eq. (D.18). We see, again, that the elements of 
Scorn provide the ratios of the output to the input state vector coefficients. Hence, Scorn is the 
scattering matrix for the combined element The quantities S™, S^, S ha , S cb , S cc , and E in 
Eq. (D.18) are, again, specified by Eqs. (D.8) and (D.14), only now the change in x a , x b , and x c 
gives elements that are different from those for the stator case. 


- 64 - 



INTERFACE 
PLANE a 


INTERFACE 
PLANE b 


INTERFACE 
PLANE b 


INTERFACE 
PLANE C 



(A) ROTOR AND ACTUATOR DISK AS SEPARATE ELEMENTS 
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(B) ROTOR AND ACTUATOR DISK AS ONE ELEMENT 


FIGURE D.2 COMBINING THE ROTOR ALONE AND ACTUATOR DISK INTO ONE ELEMENT 
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